Hall noise and transverse freezing in driven vortex lattices 
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We study driven vortices lattices in superconducting thin films. Above the critical force Fc we 
find two dynamical phase transitions at Fp and Ft, which could be observed in simultaneous noise 
measurements of the longitudinal and the Hall voltage. At Fp there is a transition from plastic flow 
to smectic flow where the voltage noise is isotropic (Hall noise = longitudinal noise) and there is 
a peak in the differential resistance. At Ft there is a sharp transition to a frozen transverse solid 
where the Hall noise falls down abruptly and vortex motion is localized in the transverse direction. 
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PACS numbers; 74.60. Ge, 74.40. +k, 05.70.Fh 

The study of the collective motion of vortex lattices in 
superconductors has brought new concepts in the non- 
equilibrium statistical physics of driven disordered media 
||l|-|l8|. The prediction |l|] of a dynamical phase transi- 
tion upon increasing drive, from a fluidlike plastic flow 
regime |^-^ to a coherently moving solid p^, has moti- 
vated an outburst of recent theoretical H-R]T experimen- 
tal ||q-[ll[, and simulation |l2-18| work. The relevant 
physics of the high velocity driven phase is controlled by 
the transverse displacements (in the direction perpendic- 
ular to the driving force) H, leading to a new class of 
driven systems characterized by anisotropic spatial struc- 
tures with transverse periodicity |^-0. Recently, these 
moving anisotropic vortex structures have been observed 
experimentally by Pardo et al. [ pl|, and their different 
features have been studied in 2D ]12|-p^ and 3D [ pTJp^ 
simulations. In this letter, we show that a better insight 
on the moving phases can be obtained from studying the 
anisotropic temporal fluctuations. We flnd two dynam- 
ical phase transitions which could be observed experi- 
mentally by measuring voltage noise ||lj,|22l both in the 
longitudinal and transversal directions. In the transverse 
direction, where the system is not driven, we can study 
diffusion, from which we deflne an effective temperature 
in analogy with equilibrium physics |]J,[7|]. 

The equation of motion of a vortex in position r^ is: 



dv- 

V^ - -^ V,C/„(r,,) - ^ V,C/p(np) + F, 



(1) 



JT^i 



where 



Ir-i — r,| is the distance between vortices 



The vortices interact with a random uniform distribution 
of attractive pinning centers with Up{r) = — Ape^^'/^^ 
with ^ being the coherence length. We normalize length 
scales by ^, energy scales by Ay, and time is normalized 
by T = rjS^'^/Ay. We consider A''^ vortices and Np pinning 
centers in a rectangular box of size L^ x Ly, and the nor- 
malized vortex density is riy = Ny£^'^ /L^Ly = B^^/^q. 
Moving vortices induce a total electric field E = — v x z. 

We study the dynamical regimes in the velocity-force 
(voltage-current) characteristics at T = 0, solving Eq.(l) 
for increasing values of F = Fy. We consider a constant 
vortex density n^ — 0.12 in a box with L^/Ly — \/3/2, 
and Ny = 64, 144, 196, 256, 400, 784 (we show results for 
Ny = 400). We take a pinning strengh of Ap/Ay = 0.2 
with a density of pinning centers Up — 0.24. We use pe- 
riodic boundary conditions and the periodic long-range 
logarithmic interaction is evaluated with an exact and 
fast converging sum [ pl| . The equations are integrated 
with a time step of At — O.Olr and averages are evalu- 
ated in 32768 integration steps after 2000 iterations for 
equilibration (when the total energy reaches a stationary 
value). Each simulation is started at F = with an or- 
dered triangular lattice and slowly increasing the force in 
steps of AF — 0.05 up to values as high as F = 8. 

First we start by looking at the vortex trajectories and 
their translational order in the steady state phases as 
was done in ||l^-[l8|]. In Figure l(a-c) we show the vor- 
tex trajectories {ri{t)} for typical values of F by plotting 
all the positions of the vortices for all the time iteration 
steps. We also study the time-averaged structure factor 
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tex i and a pinning site at Tp, ?] 



is the distance between the vor- ^'(k) = (lif Eiexp[ik • r,(t)]|2), which is shown in Fig 






is the 

Bardeen-Stephen friction and F = ^^J x z is the driv- 
ing force due to an applied current J. A two-dimensional 
superconductor is realized in thin films of thickness d 
where d <ti X, which have an effective penetration depth 
A — 2X^ /d. Since A is of the order of the sample size 
(A w 200/iTO in H), the vortex- vortex interaction is loga- 
rithmic: Uy{r) = -AylTi{r/A), with Ay = $g/87rA |,p|. 



l(d-f). In Fig. 2(a) we plot the average vortex veloc- 
ity V = {Vy{t)) = {^T,^ ^): in the direction of the 
force as a function of F and its corresponding derivative 
dV/dF. {V = E/pfJo with pf the flux flow resistivity 
and Jo = cAy/d£,^Q, therefore dV/dF is proportional to 
the differential resistivity dV/dF = p^^dE/dJ). Below 
a critical force Fc « 0.25 all the vortices are pinned and 
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FIG.l 

FIG. 1. Vortex trajectories: (a) F = 0.5, (b) F = 2.5, (c) 

F = 7. Surface intensity plot of the averged structure factor 
^(k): (d)F = 0.5, (e) F = 2.5, (f) F^7. 

there is no motion, V = 0. Above Fc, we distinguish 
between three different dynamical regimes: 

(i) Plastic flow: Fc < F < Fp. At Fc vortices start to 
move in a few filamentary channels, as was also seen in 
M. A typical situation is shown in Fig. 1(a), where a 
fraction of the vortices are moving in an intrincate net- 
work of channels. As the force is increased a higher frac- 
tion of vortices is moving. In this regime, vortices can 
move in the transverse direction (perpendicular to F) 
through the tortuous structure of channels H . We see in 
Fig. 1(d) that S'(k) only has the central peak showing 
the absence of ordering in this plastic flow regime [p|~H • 

(ii) Smectic flow: Fp < F < Ft. We observe a 
peak in the differential resistance at a characteristic force 
Fp w 1.3. At F = Fp we see that all the vortices are mov- 
ing in a seemingly isotropic channel network with max- 
imum interconnectivity In other simulations the peak 
in the differential resistance was found to coincide with 
a maximum in the number of defects |15[| and with the 
onset of orientational order ||l^. Also, the value of Fp 
was taken in the experiment of Hellerqvist et al [p| as an 
indication of a dynamical phase transition. In fact, we 
find that above Fp a new dynamic regime sets in. In this 
case, as we show in Fig. 1(b), all the vortices are moving 




FIG. 2 

FIG. 2. (a) Velocity-force curve (voltage-current character- 
istics), left scale, white points. dV/dF curve (differential re- 
sistance), right scale, black points, (b) Intensity of the Bragg 
peaks. For smectic ordering, S{Gi), Ky — : black points. 
For longitudinal ordering, S{G2,3), Ky =^ : white points. 
(c) Finite size exponent ao {S{Gi) ~ Nv"'^)- 

in trajectories that are mostly parallel to the force, form- 
ing "elastic channels" . Two small Bragg peaks appear 
in the structure factor along the Ky = axis, as seen in 
Fig. 1(e), which correspond to Gi = (±27r/ao,0). This 
is consistent with the onset of "smectic" ordering P,f2l in 
the transverse direction with elastic channels separated 

— 1/2 

by a distance r^ ag = riv ■ In this regime the trans- 
verse motion consists in vortex jumps from one chan- 
nel to another, resembling "thermally" activated tran- 
sitions induced by local chaos. The rate of these tran- 
sitions decreases with increasing force, and they corre- 
spond to the permeation modes discussed in l6l. In 
Fig. 2(b) we plot the magnitude of the Bragg peaks 
at Gi, Ss — S'(G'i), corresponding to smectic ordering 
{Ky = 0), and the other neighbouring peaks at G2 = 
±27r/ao(l/2,\/3/2) and G3 = ±27r/ao(-l/2, \/3/2) , 
Si = {S{G2) + S{G3))/2), corresponding to longitudinal 
ordering {Ky 7^ 0). We see that above Fp the intensity 
of the smectic peak Sg starts to grow and Ss ^ Si, while 
below Fp the spatial structure is isotropic, Ss = Si <^ 1. 
The Bragg peak heights depend with system size as 
S{G) ^ N~°''^ , where ctg = means long-range or- 



der (LRO), < (Tg < 1 means quasi long-range order 
(QLRO) and ac — ^ means short-range order (SRO). In 
Fig. 2(c) we show the corresponding results for the smec- 
tic peak Ss = S{Gi) for sizes iV„ = 256, 400, 784. We see 
that CTG <^ 1 m this regime: there is only short range 
smectic order, thus this phase corresponds to a hquid. In 
this sense the transition at Fp is a dynamic transition in 
the flow, as found in [^ for strong disorder {up = 1). 
When approaching Ft (see below) we see a precursor of 
QLRO (CTg < 1) but with strong fluctuations. 

(in) Frozen transverse solid: F > Ft. At a new char- 
acteristic force Ft, the jumps between channels suddenly 
stop and vortex motion becomes frozen in the direction 
perpendicular to F. An example for F > Ft is shown 
in Fig. 1(c) where we see well defined elastic channels 
parallel to F. The corresponding structure factor is in 
Fig. 1(e) where new peaks appear in ^(k) in directions 
with Ky y^ 0, like G2, G3, showing that there is some 
longitudinal ordering between the channels. These extra 
peaks are smaller than the smectic peaks, and 5'(k) is 
very anisotropic. We note in Fig. 2(a) that Ft is the 
point where the noisy behavior in dV/dF ceases. A simi- 
alar criterion was used by Bhattacharya and Higgins to 
define their dynamical phase diagram g. In Fig. 2(b) 
we see that in Ft there is an increase in the longitudi- 
nal ordering 5*;, and both Ss and Si tend to saturate at 
an almost constant value for F ^ Ft- In Fig. 2(c) we 
find that there is smectic QLRO P,0,O| with a value of 
0.5 < aa < 0.7. However, we were not able to obtain 
a reliable finite size analysis for the longitudinal peaks 
since they have large fluctuations from size to size. Also, 
we find hysteresis in Sg around Ft when decreasing F. 

A better understanding of the dynamical transitions 
can be obtained from studying the temporal behavior 
of the system in both directions. It has been observed 
experimentally that the longitudinal voltage can show 
low frequency noise [ |l9|j20| . This voltage noise reaches 
a very large value above the critical current, which has 
been attributed to plastic fiow ^^, and then the noise 
decreases for large current. In addition, even when the 
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total dc transverse voltage (I4) — (- 
it can also have fluctuations and noise p2| . In fact, 
it is easy to understand that this Hall noise will be 
closely related to the wandering and wiggling of the 
plastic flow channels and to the jumps between elas- 
tic channels in the smectic phase. We have calculated 
the power spectrum of both the longitudinal voltage, 
Sy{f) = \^ f^ dt{Vy{t)~V)exp{i2TTft)\^, and the trans- 
verse voltage, Sxif) = |y/g^dil4(t)exp(i27r/t)|2. The 
low frequency noise is defined as Px^y = lim/^o Sx.yif)- 
In Fig. 3(a) we show the values of the longitudinal noise 
Py and the Hall noise Px as a function of the force {Px,y 
were approximated from the average of the 10 lowest fre- 
quencies). We see that near the critical force, the longi- 
tudinal noise is large while the Hall noise is one order 
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FIG. 3. (a) Low frequency voltage noise vs F. Longitudinal 
voltage noise, Py, black points. Hall noise, Px, white points. 
(b) Diffusion coefficient for transversal motion, D^- The inset 
shows the average transverse quadratic displacement (A'^x) in 
the frozen phase, F > Ft. (c) Diffusion exponents (x (black 
dots) and ("^ (white dots) defined from: {[S:{t) — 5(0)]^) ~ t^^ , 
ilvi't) ~ y(0)]^) ~ t^"- The inset shows a plot of Ay defined 



fromafit ([y(t)-y(0)]^ 



Ayr for F > Ft. 



of magnitude smaller. In this regime of plastic flow the 
noise is dominated by fluctuations in the direction of mo- 
tion and by channel tortuosity H, and since there are 
few channels the Hall noise is small. When the number 
of channels increases the vortices wander more in the x 
direction and the Hall noise increases. At Fp the volt- 
age noise becomes isotropic, Px = Py This is the point 
where we have seen the highest interconnection in the 
channel network. Above Fp, the onset of elastic chan- 
nels and smectic ordering reduces the longitudinal noise. 
On the other hand, the Hall noise remains large due to 
the "activated" jumps between elastic channels. At Ft 
the Hall noise falls abruptly, nearly two orders of magni- 
tude. This corresponds to a freezing transition of vortex 
motion in the transverse direction. Above Ft there are 
no more vortex jumps between elastic channels. The low 
frequency noise can be closely related to diffusive motion 
for large times. We analyze the average quadratic dis- 
placements of vortices in both directions from their center 
of mass position {Xcm{t),Ycm{t)) as a function of time. 
We define Wx{t) = i^I]J^i(0 - 5»(0)]^> and Wy{t) = 



■W:J2^[y^it) - yM? where x^it) = x,it) - X,™(t), 
yi{t) = yi{t) — Ycm{t)- We observe that the vortex motion 
for Fc < F < F-t IS diffusive in the transverse direction, 
Wxit) ~ D^t. In Fig. 3(b) we show the behavior of the 
transverse diffusion coefheient D^- It closely follows the 
behavior of the Hall noise, D^ oc Px- The transverse 
diffusion is maximum at the peak in the differential re- 
sistance, Fp, and it has an abrupt jump to zero at Ft 
indicating the transverse freezing transition. It is inter- 
esting to note that melting transtions also show a jump 
in the diffusion coefficient. Above Ft, the transverse wan- 
dering is independent of time since vortex motion is lo- 
calized in the x direction: Wx{t) « (A^x). In the inset 
of Fig. 3(b) we show (A^a:) vs. F iov F > Ft. We 
find that (A-^a;) sa 0.02a^ at Ft, consistent with a Lin- 
demann criterion for melting 0] . Since in equibrium the 
diffusion coefficient D is proportional to T, we interpret 
Dx as an indication of an "effective temperature" Te/f. 
This means that T^jf raises above the critical force, has 
a maximum at the peak in the differential resistance and 
decreases when the motion starts to order. Following the 
same idea, above Ft we can interpret (A^x) as propor- 
tional to Teff. Here we find that (A^x) decreases with 
F, consistent with a T^ff ~ l/V behavior iQj^]- 

The long-time behavior is better understood by study- 
ing the diffusion exponents Wx (t) ~ i^^ and Wy (t) ^ i^" . 
In Fig. 3(c) we show the behavior of Cx,y We see that 
ior Fc < F < Ft, Cx ^ i corresponding to normal dif- 
fusion, while ioT F > Ft, C,x ^ corresponding to the 
freezing of transverse motion. More interestingly, the 
motion in the longitudinal direction is always superdif- 
fusive, (^y > 1. Above Ft, in the frozen phase, the longi- 
tudinal fluctuations become ballistic, C,y Ki 2. Since the 
vortex positions are localized in the x direction, we see 

that Wy{t) = {{y{t) - 2/(0))2) « Ayt\ with Ay = (A^Vy) 

the dispersion of the average velocities of the elastic chan- 
nels. In the inset of Fig. 3(c) we find that Ay decreases 
exponentially with the force Ay '^ exp(— a-F). [This 
suggests a "thermal" distribution of channel velocities, 
since Teff ^ l/F, then {A'^Vy) ~ exp(-a7Te//)]. In 
the smectic flow region an exponent \ < C,y < 2 can be 
explained from assuming that vortices move most of the 
time in the elastic channels and occasionally have diffu- 
sive jumps between channels in the transverse direction. 

In conclusion, we have obtained evidence of two dy- 
namical phase transitions which can be verified exper- 
imentally by measurements of voltage noise and Hall 
noise. The first transition at Fp is the point of isotropic 
noise and maximum transverse diffusion (i.e., maximum 
effective temperature) and corresponds to the observed 
peak in the differential resistance. The second transition 
at Ft is a freezing transition in the transverse direction, 
where the transverse diffusion vanishes abruptly and the 
Hall noise drops many orders of magnitude. 
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